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ABSTRACT 

The merger of two neutron stars results often in a rapidly and differentially rotating hypermas- 
sive neutron star (HMNS). We show by numerical-relativity simulation that the magnetic-field profile 
around such HMNS is dynamically varied during its subsequent evolution, and as a result, electromag- 
netic radiation with a large luminosity ~ 0.1B 2 R 3 Q is emitted with baryon (£?, R, and Q are poloidal 
magnetic-field strength at stellar surface, stellar radius, and angular velocity of a HMNS). The pre- 
dicted luminosity of electromagnetic radiation, which is primarily emitted along the magnetic-dipole 
direction, is ~ 10 47 (£?/10 13 G) 2 (R/10 km) 3 (^/10 4 rad/s) ergs/s, that is comparable to the luminosity 
of quasars. 

Subject headings: gamma-ray burst: general — magnetohydrodynamics (MHD) — methods: numerical 
— stars: neutron 
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1. INTRODUCTION 

Coalescence of binary neutron stars (BNS) is one of the 
most promising sources for next-generation kilo-meter 
size gravitational- wave detectors such as advanced LIGO, 
advanced VIRGO, and LCGT. Statistical studies have 
suggested that the detection rate of gra vitational waves 
emitt ed by BNS will be ~ 1-100 per year (jKalogera et al.l 
120071 ). Typical signal-to- noise ratio for such detection 
will be ~ 10 or less. Thus, it will be crucial for the detec- 
tion of gravitational waves to find electromagnetic coun- 
terparts to the gravitational- wave signals. Short-hard 
gamma-ray bursts (SGRB) h ave been inferred to accom- 
pany with the BNS merger (jNaravan et al.l 119921 : IP iran 
2004). However, this hypothesis relies on many uncertain 
assumptions; e.g., high magnetic-field strength or effi- 
cient pair annihilation of neutrino-antineutrino. In this 
article, we give a conservative estimate for the strength 
of electromagnetic signals based on a numerical-relativity 
simulation and show that a strong electromagnetic signal 
will indeed accompany with the BNS merger. 

BNS evolves due to gravitational radiation reaction 
and eventually merges. After the merger s ets in , 
there are two pos sible fates (e.g., Shi bata et al.l ((2005); 
iKiuchi et all ([2009( 1): If the total mass M is larger than 
a critical mass M c , a black hole will be formed, while 
a hypermassive neutron star (HMNS) will be formed 
for M < M c . The value of M c depends strongly on 
the equation of state (EOS) of neutron stars, but the 
latest discovery of a high-mass neutr on star with mass 
1.97±O.O4M (jDemorest et al.ll2Q10h indicates that the 
EOS is stiff and M c may be larger than the typical total 
mass of the BNS - 2.7M (|Stairsll2004l ). This indicates 
that HMNS is the like ly outcome for many B NS mergers, 
at least temporarily (Hotokez aka et al.1l2011h . 

Neutron stars in nature have a strong magnetic field 
with typical field strength at the stellar surface 10 11 - 
10 13 G. One of the neutron stars in BNS often has field 
strength smaller than this typical value as ~ 10 10 -10 11 G 
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(|Lorim er 2008}), probably because of the accretion history 
of the first-formed neutron star during the formation of 
the second one. However, at least the second one is likely 
to have the typical magnetic-field strength. 

It is reasonable to believe that each neutron star in the 
inspiral phase (before the merger sets in) has an approx- 
imately dipole magnetic field as in the isolated one. Dur- 
ing the late inspiral phase and formation of a HMNS in 
the merger phase, its magnetic-field profile will be mod- 
ified due to magnetohydrodynamics (MHD) processes. 
However, in zeroth approximation, it would be safe to 
suppose that the dipole field is dominant. Due to this 
reason, we consider the evolution of a HMNS with dipole 
magnetic fields in the following. 

One of the most important properties of HMNS is that 
it is rapidly and differentially rotating ( Shibata et al.l 
2005). The numerical simulations have shown that the 
typical angular velocity at its center is £1 J> 10 4 rad/s, 
much larger than that of ordinary pulsars, while at equa- 
torial surface it is ~ 10 3 rad/fl Because of the pres- 
ence of the differential rotation, the winding of magnetic 
fields is enhanced: Toroidal magnetic-field strength Bt 
in HMNS increases linearly with time (t) in the pres- 
ence of seed poloidal (cylindrically radial) magnetic fields 
Bp (Bt increases as ~ BpQt). The increase of the 
magnetic-field strength results in the increase of mag- 
netic pressure. Because only dilute matter is present in 
the surface of HMNS, Alfven waves are likely to propa- 
gate near the rotational axis with w <i 10 km, where 

vd = yj x 2 + y 2 , transporting electromagnetic energy 
generated in HMNS along the rotational axis; tower- type 
outflow is driven. As far as the HMNS is alive and the 
rapid rotation is present, the amplification of the toroidal 
magnetic field continues via the winding effect. Then, the 
electromagnetic energy should increase approximately as 
E B ~ BpVQ as described in Meier] (jl999h where V is an 



3 We note that even for rigidly rotating neutron stars, rotating 
magnetic field lines with a high degree of differential rotation are 
produced in the vicinity of the neutron stars because of the presence 
of a light cylinder close to their equatorial surface at ~ c/fl = 
30(f2/10 4 rad/s) -1 km with c being the speed of light. 
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effective volume for which the amplification occurfl If 
the amplification efficiently occurs near the rotation axis 
with w <; 10 km, V is approximately ~ (aR) 3 where R is 
the equatorial stellar radius ~ 10 km and a is a constant 
of O(0. 1). For a typical HMNS formed after the merger 

E B ~ 10 46 £ 2 3 Rl n 4 a 3 0A ergs/s, (1) 

where B 13 = B P /10 13 G, R 6 = R/10 km, Ol = 
ft/10 4 rad/s, and ao.i = a/0.1. We suppose a rela- 
tively high magnetic-field strength because it is likely 
to be a mplified by compress ion that occurs during the 
merger (Rez zolla et al.ll2011f ). Thus, the luminosity of 
the electromagnetic radiation will be as high as that 
of quasars for a typical field strength of a progenitor 
neutron star Bp ~ 10 12 G and of a resulting HMNS 
Bp ^ 10 13 G. If a substantial fraction of this generated 
electromagnetic energy is conv erted to electromagnet ic 
radiation (as suggested, e.g., in Naka r~fc Piranl (j2011f )). 
the merger event may be detected by telescopes as the 
electromagnetic signals. 

2. NUMERICAL SIMULATION 

Motivated by the fact mentioned above, we performed 
MHD simulation for a HMNS in general relativity. As 
the initial condition, we prepare a rapidly and differ- 
ential ly rotating HMNS in axisymmetric equilibrium 
as in IShibata et all (|2006f ): We constructed a HMNS 
model with the following piecewise polytropic EOS: P = 
^Pcoid = Kip Fl for p < p nuc and P co i d = K 2 p F2 for 
P > Pnuc- Here, P and p are the pressure and rest- 
mass density. We set Ti = 1.3, T 2 = 2.75, K x = 5.16 x 
10 14 cgs, K 2 = K lP l^- T \ and p nuc = 1.8 x 10 14 g/cm 3 . 
With this EOS, the maximum gravitational mass (rest 
mass) is 2.01M Q (2.32M ) for spherical neutron stars 
and 2.27M (2.60M©) for rigidly rotating neutron stars. 
These are similar values to those in realistic stiff EOS 
(e.g., iRead et all (j2009f ) for a review). We prepare a 
HMNS with the following physical parameters: gravi- 
tational mass M = 2.65M , baryon rest mass M5 = 
2.96M , maximum density p max — 

9.0 x 10 14 g/cm 3 , an- 
gular momentum J — 0.82GM 2 /c, ratio of polar to equa- 
torial radius 0.3, central rotation period P c = 0.202 ms, 
and rotation period at the equatorial surface 5.4P C . Here, 
G is the gravitational cons tant. The rotation law is spec- 
ified in the same way as in iBaumgarte et al.l (|2000l ) with 

the differential rotation parameter A = 0.8. This HMNS 
is similar to tha t foun d in the BNS merger simulation 
of IShibata et al.l ()2QQ5h : iKiuchi et all (|2009f ) performed 
in nuclear-theory-based EOS. In the evolution, we em- 
ploy an EOS of the form P = P co \d + (^1 — l)p(^ — £coid) 
where e is the specific internal energy and e co \d is its cold 
part determined from P co id- 

4 We note t hat the luminosity of win ds by the magnetocen- 
trifugal effect (Blandford Sz Payne 1982) has the same order of 
magnitude, but with the vertical-dominant dipole fields consid- 
ered here, this effect is not dominant. The magnetic dipole ra- 
diatio n could also play an important role as in ordinary pul- 
sars ( Vietri 1996; Lipunov & Panchenko 1996; Ioka & Taniguchi 
2000) because its luminosity oc B 2 R 6 fl 4 /c 3 (Shapiro & Teukolskv 
1983) may be comparable to the luminosity of Eq. JT} for HMNS 
with RQ/c > 0.1. However, the property of electromagnetic wave 
emission (e.g., emission direction) would be different from that we 
consider in this paper. 



Poloidal magnetic fields, for which the toroidal com- 
ponent of the vector potential has the form = 
Aqvoqvo 2 I (r 2 + tu 2 ,) 3 / 2 , are superimposed on this initial 
condition. The magnetic field in the inertial frame is 
given by B % = e % i k djAk where is the completely 
anti-symmetric tensor. Here, r 2 = w 2 + z 2 , and (tz7 , Aq) 
are constants: tz7 is chosen to be 5i? e /3-20i? e /3 where 
R e (= 12.1 km) is the coordinate radius on the equato- 
rial plane. We found that the electromagnetic luminos- 
ity shown below depends weakly on this parameter. In 
the following we show the results for tz7 = 10R e /3 and 
20R e /3. Aq determines the field strength for which we 
give the maximum magnetic-field strength 5 max « 10 13 - 
4 x 10 14 G. Here, the magnetic-field strength is defined 
by B = ^/b^W where 6 M is the 4-vector of the mag- 
netic field in the frame comoving with the fluid. With 
such strength, the magnetic pressure in the HMNS is 
much smaller than the matter pressure at its center, and 
thus, the density profile of the HMNS (except for its 
surface) is not significantly modified by the magnetic- 
field effect. Because of the presence of differential ro- 
tation, the magnetic fields may be amplified in the ex- 
ponential manner due to the magn etorotational instabil- 
ity (MRI) (jBalbus fc Hawlevlll998T ) . However, the wave- 
length for the fastest growing mode is very short (~ 10 2 - 
10 4 cm) in the present setting and it is not possible to 
resolve in numerical simulation. We here do not pay at- 
tention to the MRI but only to the winding effect. In the 
presence of the MRI effects, the electromagnetic energy 
will be increased more rapidly and the luminosity of elec- 
tromagnetic radiation may be even enhanced. Thus, this 
work would determine the lower bound of the magnetic 
luminosity, that is, however, quite high. 

MHD simulation is performed assuming that ideal 
MHD condition holds. A conservative shock captur- 
in g scheme is employed for so lving MHD equations as 
in IShibata fc Sekiguchi (2005): In the present work, a 
numerical scheme with third-order accuracy in space 
and fourth-order accuracy in time is employed. Ein- 
stein's evolution equations are solved in the fourth- 
order accuracy in space and time in the so-called 
BSSN-puncture formulation ([Shibata & Nakamura 19951 : 
Baumga rte fc Shapiro! 119 99; Campan elli et al.ll2006[ ). 

Any conservative scheme in MHD cannot handle vac- 
uum, and hence, we have to add an atmosphere of small 
density outside the HMNS. Because the matter outside 
a HMNS formed in a real merger would be dilute, the 
density of the atmosphere should be as small as possible 
to exclude spurious effects by it. We set the density of 
the atmosphere as 

J /atPmax f ^ 2_R e , /»x 

Pat ~ X /at/W(r/2i? e )-" r > 2R e , W 

where we choose n — 2 or 2.5. / at is constant, for which 
we typically give 10 -9 . We changed the values of f at from 
10- 10 to 10- 7 for £ max = 4.2 x 10 13 G and found that 
as far as B 2 / (4tt p at c 2 ) <> 1, our code works well. For a 
large value of / at , the evolution of magnetic fields is sub- 
stantially affected by the inertia of the matter. However, 
with decreasing the value of / at to J B 2 /(47rp at c 2 ) ~ 1, 
the dependence of magnetic-field evolution on the atmo- 
sphere density becomes weak, and hence, the effect of the 
artificial atmosphere does not play a role. The velocity of 
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Figure 1. Snapshots of density profiles (upper panels) and the magnetic-field strength profiles (lower panels) outside HMNS at t = 0, 2.8, 
6.2, and 12.7 ms for a model with B max = 1.7 x 10 14 G, vjq = 10R e /3, and n = 2. The region with p > 10 7 g/cm 3 are shown to be white 
in the upper panels. 



the atmosphere is set to be zero initially. With this treat- 
ment, the magnetic field is initially modified for t <i P c . 
However, such modification plays a minor role after the 
winding effect becomes dominant for the magnetic field 
amplification. We always perform simulations for a time 
much longer than P c , and focus on the stage for which a 
quasisteady state is achieved. Thus, the artificial effect 
associated with the initial setting does not matter. 

Axisymmetric numerical simulation is performed in 
the cylindrical coordinates for MHD and in the Carte- 
sian coordinates for Einstein's equation part (using the 
so -called Cartoon method). The details are described 
in lShibata fc Sekiguchil (j2005[ ) and its references to which 
th e reader may refer. Nonuniform grid is prepared as 
in iKiuchi et al I ((2008). We here impose the axial sym- 
metry to guarantee a sufficiently high grid resolution al- 
though we can perform a nonaxisymmetric simulation. 
In the present setting, the equatorial coordinate radius of 
the HMNS, R ei is covered by 150 uniform grids. Smaller 
grids with 100 and 120 were also adopted to check the 
convergence of the numerical results. Outer boundaries 
along x and z axes are located at « 170 R e « 2000 km. 

Figure [T] plots the evolution of the density profiles and 
magnetic-field strength outside the HMNS. Numerical 
simulation shows that the system evolves in the following 
manner. Because of the presence of differential rotation, 
the winding of magnetic fields proceeds. Toroidal field 
strength is increased linearly with time, and the growth 
rate is high in particular near the rotational axis (w <> 
10 km). Hence, Alfven waves propagate primarily to- 
ward ^-direction near the rotational axis, and magnetic- 
field strength there also increases. After the substantial 
winding, the magnetic pressure B 2 /8tt becomes larger 
than the gravitational potential energy density near the 
polar surface, ~ GMp/H where H(= 0.3R e ) is the verti- 
cal coordinate radius of the HMNS. Then, the matter of 
the HMNS in the vicinity of its polar surface is stripped, 
and an outflow is driven. Because the magnetic energy 
density B 2 /4ir is comparable to or slightly smaller than 



the rest-mass density - (H/GM)(B 2 /8tt) > B 2 /4ttc 2 , 
the outflow is mildly relativistic with the velocity of or- 
der 0.1c in the vicinity of the HMNS. However, in the 
region far from the HMNS, the outflow velocity could 
be fairly relativistic ~ 0.9c near the rotation axis (see 

Fig. E]). 

Figure [2] plots evolution of the matter and electromag- 
netic energy ejection rates Em and Eb defined by 



E 



(3) 



where we substitute the stress energy tensor for the mat- 
ter and electromagnetic fields into T r £ , respectively, g 
is the determinant of the spacetime metric. The surface 
integral is performed for r ~ 480 km. We checked that 
the luminosity depends only weakly on the radius of the 
surface integral. The top-left and top-right panels show 
the evolution of Em and E B for £? max = 4.2 x 10 13 G, 
and the bottom- left panel show the evolution of Eb for 
5 max = 1.7 x 10 13 -4.2 x 10 14 G. The bottom-right panel 
shows the ejection rates integrated only for the angle 
AO = tt/20 for £ max = 4.2 x 10 13 G. The top-left panel 
is for wo = 10i? e /3, and n = 2 with three grid resolu- 
tions, and the top-right panel is for different values of 
wo and n. The top- left and top-right panels show that 
irrespective of the grid resolution, wq : and n, the total 
ejection rates are 

E M ~10 48 B 2 3 R 3 6 n 4 ergs/s, (4) 

E B ~ 10 47 £? 3 i^fi 4 ergs/s. (5) 

The bottom-left panel indeed shows that the scaling re- 
lation with respect to the magnetic-field strength holds 
(the same scaling also holds for Em)- Another important 
point is that these ejection rates do not significantly vary 
in time. Thus, a quasisteady outflow is driven. 

The value of Eb agrees approximately with the pre- 
diction of Eq. (pQ), implying that the scenario described 
in Sec. 1 is correct. The ratio of Em/Eb is of order 
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Figure 2. Top-left: Total matter energy (upper curves) and magnetic energy (lower curves) ejection rates as functions of time for 
B max = 4.2 x 10 13 G, vuo = 10Re/3, and n = 2. The results with three different grid resolutions are plotted: The solid, dashed, and dotted 
curves show the results with high, middle, and low resolutions. Top-right: The same as the top-left panel but for (vjQ,n) = (10-R e /3, 2) 
(solid curves), (KXR e /3, 2.5) (dashed curves), and (20-R e /3,2) (dotted curves). Bottom-left: Evolution of E B for B max = 1.7 x 10 13 , 
4.2 x 10 13 , 1.7 x 10 14 , and 4.2 x 10 14 G. Bottom-right: The same as the top panels but for the amount integrated only for the angle 
AO = 7r/20 radian from the rotation axis for (mo,n) = (10i? e /3, 2) (solid curves) and (20i? e /3, 2) (dotted curves). The thick curves denote 
the electromagnetic luminosity. For all the figures, r denotes the radius for which the surface integral of Eq. ((3} is carried out. 



2c 2 H/GM ~ 10. This agrees approximately with the 
value required for the mass stripping. 

Comparison among top- left, top-right, and bottom- 
right panels of Fig. [2] shows that the electromagnetic 
energy is mainly emitted in the direction near the ro- 
tation axis. By contrast, the matter energy is emitted 
in a fairly isotropic manner. Along the rotational axis, 
Em/Eb is of order unity (~ 1-10). This is the reason 
that the outflow along the rotation axis can be mildly 
relativistic. 

The amount of angular momentum loss by the mat- 
ter ejection and electromagnetic radiation in the time 
duration At = P c is much smaller than the total angu- 
lar momentum because of our choice of magnetic-field 
strength <J 3 x 10 14 G. This implies that the matter 
and electromagnetic waves are continuously ejected, and 
a quasisteady outflow is formed. Figure [3] plots the 
density and velocity profiles along the rotation axis for 
£ max = 4.2 x 10 13 G, w = 10i? e /3, and n = 2. The den- 
sity profile (left panel) indeed shows that the averaged 
density does not change significantly with time. The den- 
sity decreases with the radius. In these examples, the 
power-law index is roughly n p ~ 1.5-2 (p oc r _np ), but 
this number depends on the initial setting for £? max and 
wq and varies with time. The velocity profile (right panel 



of Fig. [3]) shows that the outflow is mildly relativistic. 
The profile varies in a short time scale. The maximum 
velocity is ~ 0.9c as mentioned above: This maximum 
depends weakly on the setting of the atmosphere; for the 
lower atmosphere density, the maximum speed is larger. 
The averaged magnitude of the outflow velocity in time 
is ~ 0.4-0. 6c (which also depends weakly on the setting 
of the atmosphere). Because the corresponding Lorentz 
factor of the jet is Tj < 2, the relativistic beaming effect 
(i.e., observable viewing angle is < l/Tj) is relatively 
small, and thus, this jet may be observable from a large 
solid angle. 

3. DISCUSSION 

The centrifugal force due to rapid and differential ro- 
tation plays a crucial role for supporting strong self- 
gravity of HMNS. This suggests that any HMNS will 
eventually collapse to a black hole after a substantial 
loss of angular momentu m by gravitational-wave emis- 
sion (Shi bata et al.l 120051 ) and/or after a substantial an- 
gular momentu m transport insid e it due to magnetic 
viscous effects (jDuez et al.l [2006) in a realistic situa- 
tion. The predicted lifetime is of order 10-100 ms. 
After the collapse to a black hole, a system com- 
posed of a black hole and compact accretion disk 
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Figure 3. Snapshots of density (left panel) and velocity profiles (right panel) along the rotation axis for £> max = 4.2 x 10 13 G, vjq = 10R e /3 : 
and n = 2. At t = 0, the density is determined by the atmosphere and decreases in the form p oc r -2 for the large radius. 



will be formed dShibata et all [20061: iDuez et all 120061 : 
iRezzolla etHll2QlH iHotokez akaet alJl2011h ."Then. the 
electromagnetic radiation is likely t o be emitted through 
the B l andford-Znajek mecha n ism (Blandford & Z naiekl 
[19771 : iMcKinnev fc Gammid 120041 : iMcKinnevI 12005). 
However, the lifetime of the accretion disk will not be 
longer than ~ 100 ms because of the viscous evolution. 
Hence, the magnetic energy generation and mass ejection 
are likely to continue only for a short time scale of order 
100 ms after the BNS merger. 

However, the luminosity is quite high. If the kinetic 
energy of the matter and/or electromagnetic energy are 
efficiently converted to an electromagnetic signal, it will 
be observed even if the merger happens at a distance of 
several hundred Mpc as discussed, e.g., in Naka r~fc Piranl 
(|2011[ ). Strong radio afterglow emission could be ex- 
pected when the jet propagates in the matter that may 
be ejected from neutron stars and HMNS during the 
merger. The ejected mass that can be estimated as 
E M AT/c 2 ~ 1O- 6 M ^ 3J R^ 4 ( AT / O .1 s) could yield 
radioa ctive elements and b e observed like dim super- 
novae (jLi fc Paczv hski 1998|). 

In this article, we focus only on the HMNS with con- 
servative magnetic field strength B ~ 10 13 G (achieved 
by compression of ordinary field strength ~ 10 12 G). If 
the m agnetic field strength were a s high as that of mag- 
net ar (jWoods fc Tho mpson 2006), i.e., B ~ 10 15 G, the 
electromagnetic luminosity would reach 10 51 ergs/s (see 
the bottom- left panel of Fig. [2j). This value is as high as 
the luminosity of GRBs, and because the expected time 
duratio n is less than 1 s, this is also a candidate model for 
SGRB (Nakar 2007). The canonical peak isotropic lumi- 
nosity of SGRB is rsj 10 51 ergs/s, which is consistent with 
this estimation, assuming the jet solid angle as Qj ~ 0.1 
and the conversion efficiency from the jet kinetic energy 
into the gamma rays as r] ~ 0.1. Hence, if one of neutron 
stars in BNS has a large magnetic-field strength or the 
magnetic-field strength is significantly amplified during 
the merger process or in the formed HMNS, this may be 
observed as a SGRB. 
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